Hadron production by quark combination in central Pb+Pb collisions at ^/s]^ =17.3 GeV 
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The quark combination mechanism of QGP hadronization is applied to nucleus-nucleus collisions at top SPS 
energy. The yields, rapidity and transverse momentum distributions of identified hadrons in most central Pb-l-Pb 
collisions at ^snn = 17.3 GeV are systematically studied. The calculated results are in agreement with the 
experimental data from NA49 Collaboration. The longitudinal and transverse collective flows and strangeness 
of the hot and dense quark matter produced in nucleus-nucleus collisions at top SPS energy are investigated. It 
is found that the collective flow of strange quarks is stronger than light quarks, which is compatible with that at 
RHIC energies, and the strangeness is almost the same as those at -\fsNN = 62 A, 130, 200 GeV. 
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I. INTRODUCTION 



Lattice QCD predicts that at extremely high temperature 
and density, the confined hadronic matter will undergo a 
phase transition to a new state of matter called quark gluon 
plasma (QGP) llIlQl- The relativistic heavy ion collisions 
can provide the condition to create this deconfined partonic 
matter In general, two approaches are used to study 
the properties of the deconfined hot and dense quark mat- 
ter produced in AA collisions. One is studying the high pr 
hadrons from initial hard jets, in which one can recur to 
the perturbative QCD to a certain degree |@]. The other is 
investigating the properties of thermal hadrons frozen out 
from the hot and dense quark matter. For the latter, the 
hadronization of the hot and dense quark matter (a typi- 
cal non-perturbative process) is of great significance. Only 
through a reliable hadronization mechanism, can we re- 
versely obtain various information of QGP properties from 
the final state hadrons measured experimentally. The abun- 
dant experimental data jsj 0] and phenomenological stud- 
ies 10, i, i, [H El at RHIC energies suggest that 
quark combination mechanism is one of the most hopeful 
candidates. The two most noticeable results are the suc- 
cessful explanation of the high baryon/meson ratios and 
the constituent quark number scaling of the hadronic el- 
liptic flow in the intermediate transverse momentum range 
lH, 1^, which can not be understood at all in the partonic 
fragmentation picture. Recently, the NA49 Collaboration 
have measured the elliptic flow of identified hadrons at top 
SPS energy ifisll . and found that the quark number scal- 
ing of elliptic flow was shown to hold also. It immedi- 
ately gives us an inspiration of the applicability for the 
quark combination at top SPS energy. On the other hand, 
the NA49 collaboration have found three interesting phe- 
nomena around 30A GeV ifTill . i.e. the steepening of the 
energy dependence for pion multiplicity, a maximum in 
the energy dependence of strangeness to pion ratio and a 
characteristic plateau of the effective temperature for kaon 
production. These phenomena are indicative of the onset 
of the deconfinement at low SPS energies. One can es- 
timate via Bjorken method that the primordial spatial en- 
ergy density in Pb+Pb collisions at top SPS energy is about 



3.0 GeV//OT^ ifTsll . exceeding the critical energy density 
(about 1 GeV/fm^) predicted by Lattice QCD. Therefore, 
the deconfined hot and dense quark matter has been prob- 
ably created, and we can extend the quark combination 
mechanism to SPS energies. 

As is well known, hadron yield is one of the most ba- 
sic and important observables which can help us to test the 
understanding of the hadronization mechanism for the hot 
and dense quark matter created in the relativistic heavy ion 
collisions. In most of recombination/coalescence models, 
the hadron wave function is necessary to get the hadron 
yield. As the wave functions for almost all hadrons are un- 
known at present, it is difficult for these models to study 
this issue quantitatively In addition, these mod- 

els do not satisfy the unitarity which is important to the 
issue as well jlTll . Different from those models, the quark 
combination model ifisi [l9tl uses the near-correlation in 
phase space and SU/(3) symmetry, instead of hadron wave 
function, to determine the hadron multiplicity. In addi- 
tion, the model satisfies unitarity as well and has repro- 
duced many experimental data at RHIC ll20l l2ll l22l l23ll . 
Therefore, we apply it in this paper to systematically study 
the yields, rapidity and transverse momentum distributions 
of various hadrons in most central Pb-i-Pb collisions at 
yfsNN - 17.3 GeV. On one hand, one tests the applica- 
bility of the quark combination mechanism at this colli- 
sion energy. We note that the first attempt of the mech- 
anism at SPS energies, for hadron yields alone, was the 
ALCOR model |24, 25]. Now, the rich experimental data 
of hadron multiplicities and momentum spectra provide an 
opportunity to make a further, systematical and even deci- 
sive test of the mechanism at SPS energies. On the other 
hand, the parton momentum distributions at hadronization, 
which carry the information on the evolution of the hot and 
dense quark mater, are extracted from the final hadrons at 
top SPS energy and compared with those at RHIC ener- 
gies. We concentrate the comparison on two properties 
related to strange hadron production. One is the difi^er- 
ence in collective flow between light and strange quarks, 
which occurs at RHIC energies 1231 1261 . The other is the 
strangeness enhancement, a significant property of QGP 
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The paper is arranged as follows. In the next section, we 
make a brief introduction to the quark combination model. 
In section III, we calculate the yields and rapidity distri- 
butions of identified hadrons in most central Pb+Pb col- 
lisions at -s/snn - 17.3 GeV. In section IV, the results 
of transverse momentum distributions of various hadrons 
are shown. In section V, we firstly make a detailed analy- 
sis of the longitudinal and transverse collective flow of the 
hot and dense quark matter at top SPS energy. Secondly, 
the energy dependence of the strangeness in the hot and 
dense quark matter is extracted and analyzed. Section VI 
presents summary. 



II. AN INTRODUCTION OF THE QUARK 
COMBINATION MODEL 



The starting point of the model is a color singlet sys- 
tem which consists of constituent quarks and antiquarks. 
All kinds of hadronization models demand that they sat- 
isfy rapidity or momentum correlation for quarks in the 
neighborhood of phase space. The essence of this corre- 
lation is in agreement with the fundamental requirement 
of QCD ll28ll . According to QCD, a qq may be in a color 



octet or a singlet. The color factors {(qq)s 



6 

and {(qq)i\-^^j^\(qq)i) - which means a repulsive or 
an attractive interaction between them. Here A" are the 
Gell-Mann matrices. If they are close with each other in 
phase space, they can interact with sufficiently time to be 
in the color singlet and form a meson. Similarly, a qq 
can be in a sextet or an anti-triplet, and the color factors 
<(W)6|^I(W)6> = 5 ^nd {(qq)i\^\(qqh) = If its 
nearest neighbor in phase space is a q, they form a baryon. 
If the neighbor is a q, because the attraction strength of the 
singlet is two times that of the anti-triplet, qq will win the 
competition to form a meson and leave a q alone to com- 
bine with other quarks or antiquarks. Based on the above 
QCD and near-correlation in phase space requirements, 
we had proposed a quark combination rule(QCR) ifisl l28ll 
which combines all these quarks and antiquarks into initial 
hadrons. When the transverse momentum of quarks are 
negligible, all q and q can always line up stochastically in 
rapidity. The QCR reads as follows: 

1 . Starting from the first parton (q or q) in the line; 

2. If the baryon number of the second in the line is of 
the different type from the first, i.e. the first two par- 
tons are either qq or qq, they combine into a meson 
and are removed from the line, go to point 1 ; Other- 
wise they are either qq or qq, go to the next point; 

3. Look at the third, if it is of the dififerent type from the 
first, the first and third partons form a meson and are 
removed from the line, go to point 1 ; Otherwise the 
first three partons combine into a baryon or an anti- 
baryon and are removed from the line, go to point 
1. 



The following example shows how the above QCR 
works 

-> M{qiq2) B(q-^q^q^) M(q(,qj) B(q^qgqio) Miq^^qn) 

M{q\{qxs) B{qi4qi(,qn) ^(^18^19^20) 

If the quarks and antiquarks are stochastically arranged 
in rapidity, the probability distribution for pairs of 
quarks and antiquarks to combine into M mesons, B 
baryons and B anti-baryons is 



^^,^_ 2NimnM + 2B-iy. ^^_,^ 

^^^^^^ (2N)lMl(Blf ^ 



(1) 



Hadronization is the soft process of the strong interac- 
tion and is independent of flavor, so the net flavor number 
remains constant during the process. In the quark com- 
bination scheme, this means that the quark number for 
each certain flavor prior to hadronization equals to that 
of all initially produced hadrons after it. Obviously the 
quark number conservation is automatically satisfied in the 
model. It is different from the non-linear algebraic method 
in ALCOR model ll24ll where normalization factor for each 
quark flavor is introduced with the constraint of the quark 
number conservation. 

The average number of initially produced mesons M{N) 
and baryons B(N) are given by 



(MiN)) = J]J]mXmb(N) 

M B 

{B(N)) = J]J]^Xmb(N). 



(2) 
(3) 



Then the multiplicity of various initial hadrons is obtained 
according to their production weights 



<m;"""") = Cm,{M(N)), {BY""') = Cb,{B{N)), (4) 

where Cm, and Cb, are normalized production weights for 
the meson Mj and baryon Bj, respectively. If three quark 
flavors are considered only, we can obtain the production 
weights using the SU f (3) symmetry with a strangeness 
suppression factor As Usi which are listed in Table 
H] The extension of the symmetry to excited states, ex- 
otic states and more quark flavors is also straightforward 
SHI HI- 

Considering the decay contributions from the reso- 
nances, we can obtain the yields of final state hadrons 



initial 



> + y B,{j ^ i){hj}. 



J 



(5) 



where the B, ( /' i) is the weighted decay branching ratio 
for hj to hi 1 30]. 

In principle, the hadron production probability should 
be calculated from the matrix element {qq\M} for me- 
son or {qqq\B) for baryon. However, the wave functions 
for almost all hadrons which are governed by the non- 
perturbative QCD are unknown at present. It is difficult 
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TABLE I: The normalized production weight for baryons and 
mesons in the SU/-(3) ground state, r, is the number of strange 
quarlcs in hadron. The ratio of the vector (/'' = 1") to pseu- 
doscalar (/'' =0") meson follows the spin counting, while that 
of the decouple! {J'' = | ) to octet {J'' = 5 ) baryon suffers a 
spin suppression effect; see Ref. fT9l . l29h for details. 
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except 
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Cb, 
Ca 



(2+,i.,)3(2y,+i) 
CjP = Q.o 



/ij , except 

= Ca(1520) = 



2(2+.i,-)-' 



to Study the production of hadrons quantitatively through 
their wave functions. In view of this, the hadron pro- 
duction probability in our model is determined by the 
SU/-(3) symmetry with a strangeness suppression. This 
symmetry has been supported by many experiments, par- 
ticularly by the coincidence of the observed As obtained 
from various mesons and baryons ll3lll . Therefore, the 
model can quantitatively describe many global properties 
for the bulk system by virtue of the Monte Carlo method 

||il|23,|2ll|2i|2ll2i. 

When applying the model to describe the hadronizaton 
of the hot and dense quark matter produced in heavy ion 
collisions, the net-baryon quantum number of the system 
perplexes the analysis formula of Eq. [T]but it can be eas- 
ily evaluated in Monte Carlo program. On the other hand, 
the the transverse momentum of quarks is not negligible 
due to the strong collective flow of quark matter. In princi- 
ple, we should define the QCR in three-dimensional phase 
space, but it is quite complicated to have it because one 
does not have an order or one has to define an order in 
a sophisticated way so that all quarks can combine into 
hadrons in a particular sequence. In practice, the combi- 
nation is still put in rapidity and meanwhile the maximum 
transverse momentum diff'erence between (anti)quarks 
are constrained as they combine into hadrons. The trans- 
verse spectra of hadrons have a relationship with the quark 
spectra as follows (e.g. for meson) 



cINm 



""J '^^Pi.T'^^P2,t/9(Pi,t)/?(P2,t)'^^(Pt - Pl,T ' Plj) 



X &(Ap 



■ P2.tI) 



(6) 



where /^/^(Pj) is the transverse momentum distribution of 
the quark/antiquark, assumed to be rapidity-independent in 
present work. The superscript asterisk denotes the quark 
momentum in the center-of-mass frame of formed hadron. 
The limitation Ap is treated as parameter in our study, and 
fixed to be A^ — 0.3 GeV for mesons and Ap = 0.6 GeV 
for baryons both at RHIC and SPS energies. Note that the 
spectrum normalization is determined by the multiplicity 
in Eq. [H i.e. the constraint of the parameter Ap on the 
hadron yield is neglected. 

One issue that is often questioned is the energy and en- 
tropy conservation in quark combination process. As the 
non-perturbative QCD is unsolved, there is no rigorous 



theory which can incorporate the partonic phase as well as 
hadronic phase, thus it is difficult to justify or condemn this 
issue in essence at the moment. As we know, a lot of the 
experimental phenomena in intermediate transverse mo- 
mentum range at RHIC can be explained beautifully only 
in the quark combination scenario. It suggests that maybe 
this 'puzzling' issue does not exist. As far as the quark 
combination itself is concerned, there is no difference for 
the combination occurred in the different (intermediate or 
low) transverse momentum range. Therefore, whether the 
properties of low pj hadrons can be reproduced or not is 
also a significant test of the quark combination mechanism, 
as the vast majority of hadrons observed experimentally 
are just these with low transverse momentum. 



III. HADRON YIELDS AND RAPIDITY 
DISTRIBUTIONS 

In high energy nucleus-nucleus collisions, the energy 
deposited in the collision region excites large numbers of 
newborn quarks and antiquarks from the vacuum. Sub- 
sequently, the hot and dense quark matter mainly com- 
posed of these newborn quarks will expend hydrodynami- 
cally until hadronization. The net-quarks from the collid- 
ing nuclei still carry a fraction of beam energy, thus their 
evolution is different from the newborn quarks. One part 
of net-quarks are stopped in the hot and dense quark mat- 
ter, and hadronize together with it. The other part of net- 
quarks penetrate the hot quark matter, and run up to the 
forward rapidity region. The latter, together with small 
amount newborn quarks, form the leading fireball. Their 
hadronization should be earlier than that of the hot and 
dense quark matter with a prolonged expansion stage, and 
the hadronization outcomes consist of nucleons and small 
mount of mesons. 

The current version of quark combination model sim- 
ulates only the hadronization of the hot and dense quark 
matter and subsequently decays of resonances. One in- 
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FIG. 1: (Color online) Rapidity spectra of newborn quarks and 
net-quarks at hadronization in most central Pb+Pb collisions at 
= 17.3 GeV. 
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dispensable input is the momentum distributions of ther- 
mal quarks and antiquarks at hadronization, which are the 
results of the hydrodynamic evolution in partonic phase. 
In order to focus attentions on the test of quark combi- 
nation mechanism in this and next sections, we reversely 
extract the quark distributions by fitting the experimental 
data in the model. A detailed analysis of quark distribu- 
tions at hadronization will be in section V. The solid and 
dashed lines in Fig. [T] show the rapidity distributions of 
newborn light and strange quarks at hadronization respec- 
tively, obtained from the and data ll32ll .The dotted- 
dashed line is the rapidity distribution of net-quarks in the 
hot and dense quark matter, which is extracted from net- 
proton data ifssll . 

Firstly, we calculate the yields and rapidity densities at 
midrapidity of various hadrons in most central Pb-i-Pb col- 
lisions at yfsNN - 17.3 GeV. The results are shown in Ta- 
ble From the energy dependence of the rapidity density 
for net-baryon [i39[l , one can see that the nucleus-nucleus 
collisions at SPS energies exhibit a strong stopping power 
Therefore, the leading particles contribute little to yields 
and rapidity distributions of various hadrons. 

The calculated yields and rapidity densities of vector 
meson are shown to be about twice as high as the exper- 
imental data. The results of other hadrons are basically in 
agreement with the experimental data, but slight deviations 
exist also. The overpredictions of ^ meson may be asso- 
ciated with the exotic particle /o(980), which has a possi- 
ble tetraquark structure containing a strange quark and a 
strange antiquark ll40ll . As a bond state containing strange 
components, it has a slightly lower mass than (p meson but 
is not included in the SU/(3) ground states. In the present 
work, we consider only the production of 36 - plets of me- 
son and 56 - plets of baryon in the SU/(3) ground states, 
and the excited states and exotic states are not taken into 
account. The /o(980) multiplicity is found to be nearly the 
same as meson in the e^e" annihilations | (3Q[| . The nij 
distribution of /o(980) measured by STAR Collaboration 
in Au+Au collisions at y^wy - 200 GeV is also shown to 
be comparable to that of (p IMtI |42|i . Therefore, the over- 
prediction of (p meson can be removed by incorporating the 
/o(980) production. 

Subsequently, we will calculate the longitudinal rapidity 
distributions of various hadrons. Due to the deviations in 
hadron yields, it is difficult to directly compare the calcu- 
lated hadron spectra with the experimental data. In order 
to focus attentions on the property of hadron momentum 
spectra, we will scale the calculated rapidity densities to 
the center value of the experimental data when we show 
the hadronic rapidity and pT spectra in Fig. |2]and|4]respec- 
tively, thereby removing these deviations in hadron yields. 

Pion is the lightest and most abundant hadron produced 
in AA collisions, and its momentum distribution can best 
reflect the global evolution property of the hot and dense 
quark matter In various models of high energy heavy ion 
collisions, the reproduction of pion meson is always taken 
as a paramount test of models. In Landau hydrodynamic 
model [14311 . the rapidity distribution of pion can be well 



TABLE II: The yields (left) and rapidity densities at midrapid- 
ity (right) of identified hadrons in central Pb+Pb collisions at 
■JJnn = 17.3 GeV. The experimental data are taken from Ref. 
jslii, 35, 36, 37, 38] 
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described and the sound velocity (which is an important 
physical quantity standing for the property of the hot and 
dense quark matter) can be extracted from the pion distri- 
bution. For other hadrons, such as kaons, protons. A, H 
and so on, the Landau model can not describe their rapid- 
ity distributions with the same sound velocity or freeze- 
out temperature jiii l45l l46ll . For a systematic descrip- 
tion of the rapidity distributions of various hadrons, the 
detailed longitudinal dynamics, e.g. the evolution of net- 
baryon density which will result in the yield and spectrum 
asymmetry between hadron and antihadron, should be in- 
cluded. In addition, the hadronization mechanism is espe- 
cially important to describe the differences in the yield and 
momentum distribution of various hadron species. Using 
the extracted quark distributions in Fig. [T] we calculate the 
rapidity distributions of pions, kaons, A(A), H"(H*) and 
(p in most central Pb+Pb collisions at yfstm - 17.3 GeV. 
The results are shown in Fig. |2]and are compared with the 
experimental data. The calculated rapidity spectrum of (f) 
meson is narrow than the latest data of NA49 Collaboration 
(open circles in the second panel), but is in good agreement 
with previous data (filled circles). The rapidity spectra of 
other hadrons are well reproduced. One can see that the 
quark combination mechanism is applicable for describing 
the longitudinal distributions of various hadrons at top SPS 
energy. 



IV. HADRON TRANSVERSE MOMENTUM 
DISTRIBUTIONS 

In this section, we calculate the transverse momentum 
distributions of various hadrons in the midrapidity range. 
In this paper, we only consider the hadronization of the hot 
and dense quark matter. The transverse momentum invari- 
ant distribution of constituent quarks at hadronization is 
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FIG. 2: (Color online) The scaled rapidity distributions of identified hadrons in most central Pb+Pb collisions at -sJsnn = 17.3 GeV. 
The contributions from leading particles are not included. The open circles of ip data in the second panel are the latest results measured 
by NA49 Collaboration jssll . and filled circles are the previous ones iBtIi . Other experimental data are taken from Ref. |[32l [37l[3^ . 
The open symbols of A^°, A and H show data points reflected around midrapidity. 



taken to be an exponential form expi-jnr/T), where T is 
the slope parameter which is also called effective produc- 
tion temperature. Fig. [3] shows the midrapidity pj- spec- 
tra of constituent quarks at hadronization in most central 
Pb-i-Pb collisions at y/s^N - 17.3 GeV. The spectra of 
newborn light and strange quarks are extracted from the 
data of 7T^ and respectively ll48ll . The net-quark distri- 
bution is fixed by the data of /K^ ratio as a function of 
Pt iH]. A detailed analysis of the quark pr spectra will 
be shown in section V. 




12 3 

PtLGbV/c] 



FIG. 3: (Color online) The transverse momentum distributions 
of constituent quarks in the midrapidity region at hadronization 
in most central Pb+Pb collisions at y5jviv = 17.3 GeV. 

Fig. m shows the calculated pr spectra of pions, kaons, 
protons, A(A), S"(S*) and <f> in most central Pb+Pb col- 
lisions at yfsNN - 17.3 GeV. For kaons, protons, A and 
f2, the spectral slopes of antihadrons measured experimen- 
tally are all steeper than those of hadrons iHl Hi, Isoll . 
However, the spectrum of is abnormally steeper than 
that of E* (yf]. Our predictions are in good agreement 
with all the data except . 

The exponential function e\p{-mx/T) is often used ex- 
perimentally to fit the transverse momentum distributions 
of identified hadrons in the low pr range, and to extract 



the effective production temperature T of various hadrons. 
It is found at top SPS energy that all final-state hadrons 
except pion meson have much higher T than the critical 
temperature ifsill . which indicates a strong collective flow 
at this collision energy. It is regarded in Ref. lf52ll that 
this observed flow mainly develops in the late hadronic 
rescattering stage. But results in Fig. |2] and Fig. |4] all 
show that both longitudinal and transverse spectra of vari- 
ous hadrons can be coherently explained by the same quark 
distributions, respectively. It suggests that the observed 
flow should mainly come from the expansive evolution 
stage of the hot and dense quark matter before hadroniza- 
tion, but not from the post-hadronization stage. In addi- 
tion, the same constituent quark spectra contained in light, 
single- and multi- strange hadrons also imply that the hot 
quark matter hadronize into these initial hadrons almost at 
the same time, i.e. the hadronization is a rapid process. 

In Fig. |5] we calculate the pr spectra of pions, kaons, 
protons, A(A), S"(H*), <f) and f2(Q) in most central Au+Au 
collisions at top RHIC energy. The momentum distribu- 
tions of constituent quarks at hadronization are taken to 
be exp(-m7/0.375) for strange quarks and exp(-m7-/0.34) 
for light quarks. The numbers and rapidity spectra of the 
light and strange quarks and antiquarks have been obtained 
in the study of the longitudinal hadron production ifssll . 
Here, the result of (p meson is multiplied by a factor 0.5. 
One can see that the pr spectra of various hadrons are well 
reproduced. 

The baryon/meson ratio as a function of pr is sen- 
sitive to the hadronization mechanism. As we know, 
the observed high baryon/meson ratios in the interme- 
diate Pt range at RHIC energies ijst] can not be under- 
stood at all in the scheme of parton fragmentation, but 
can be easily explained in the quark combination mech- 
anism. Fig. |6] shows the model predictions of 'p/n^, 
A//rf and at midrapidity in both central Pb+Pb colli- 
sions at -\/snn = 17.3 GeV and central Au+Au collisions 
at yfs^ = 200 GeV. In the intermediate pj range where 
the hadron production is dominated by the combination of 
thermal quarks, the baryon/meson ratios increase with the 
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FIG. 4: (Color online) The scaled transverse momentum distributions of identified hadrons at midrapidity in most central Pb+Pb 
collisions at sfsNN = 17.3 GeV. Only the combination of thermal quarks is taken into account. Solid lines are the calculated results of 
hadrons and dashed lines for antihadrons. The experimental data are from Ref. ll35Ll37Ll48ll 




2 3 

P.r[GeV/c] 



2 3 

PT-[GeV/c] 



2 3 

PT-[GeV/c] 



FIG. 5: (Color online) The transverse momentum distributions of identified hadrons at midrapidity in most central Au+Au collisions at 
-\fsNN = 200 GeV. Only the combination of thermal quarks is taken into account. Solid lines are the calculated results of hadrons and 
dashed lines for antihadrons. The experimental data are from Ref. f^ [53l[54ll 



increasing pr- One can see that the experimental data in 
this region are well reproduced. The falling tendency of 
measured baryon/meson ratios after peak position is owing 
to the abundant participation of jet quarks, which is beyond 
the concern of the present paper Besides the hadronization 
mechanism, the baryon/meson ratio in the intermediate pr 
region is also influenced by other two factors. One is the 
nuclear stopping power in collisions. Comparing with the 



strong collision transparency at top RHIC energy ||39tl, the 
strong nuclear stopping at top SPS energy causes the de- 
tention of abundant net-quarks in the midrapidity region. 
These net-quarks significantly suppress the production of 
anti-baryons and enhance that of the baryons. Therefore, 
the p/tt" ratio at top SPS energy is much lower than that at 
top RHIC energy while the A/K'^ ratio at top SPS energy 
is higher than that at top RHIC energy. The other is the 
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PT[GeV/c] PT.[GeV/c] PT[GeV/c] 

FIG. 6: (Color online) The ratios ofp/n^, A/K^ and Q./(f) at midrapidity in most central Pb+Pb collisions at y%iv = 17.3 GeV and 
Au+Au collisions at 200 GeV. Only the combination of thermal quarks is taken into account. Solid lines are the calculated results in 
Au+Au collisions and dashed lines for Pb+Pb collisions. The experimental data are from Ref. (M l48Ll54ll55l57ll 



momentum distribution of constituent quarks at hadroniza- 
tion. This can be illustrated by ratio because the pro- 
duction of these two hadron species is less influenced by 
the net-quarks. The calculated Q./(p ratio shows a weak 
dependence on the collision energy in the intermediate pr 
range. The well description of various baryon/meson ra- 
tios in such a wide energy range is an indication of the 
universality for the quark combination mechanism. 



V. ANALYSIS OF PARTON DISTRIBUTIONS AT 
HADRONIZATION 

The constituent quark distributions at hadronization 
carry the information on the evolution of the hot and dense 
quark matter in partonic phase. In this section, we focus at- 
tentions on the longitudinal and transverse collective flows 
and strangeness enhancement of the hot and dense quark 
matter produced at top SPS energy. 



A. The longitudinal and transverse collective flow 

Due to the thermal pressure, the hot and dense quark 
matter created in high energy heavy ion collisions will ex- 
pand during the evolution before hadronization. The longi- 
tudinal and transverse collective flow of final hadrons mea- 
sured experimentally is the exhibition of this early thermal 
expansion in the partonic phase. Utilizing the relativistic 
hydrodynamic evolution of the hot and dense quark matter, 
one can obtain the collective flow in quark level from the 
extracted quark momentum distributions, and compare it 
with that at RHIC energies. 

There are two well known hydrodynamic models for 
the description of the space time evolution of the hot and 
dense quark matter produced in heavy ion collisions. One 
is Bjorken model 15811 which supposes that the collision 
is transparent, and it is appropriate to extremely high en- 
ergy collisions, such as LHC. The other is Landau model 
14311 with an assumption of the full stopping for nucleus- 



nucleus collisions. The longitudinal evolution result is 
equivalent to the superposition of a set of thermal sources 
in rapidity axis, with a (Bjorken) uniform or (Landau) 
Gaussian weight. In general, when applying the model 
to describe the hadron rapidity distributions, different pa- 
rameter values are required to make a good fit of different 
hadron species liill . In this paper, we apply the hydro- 
dynamic description to the evolution in quark level, thus 
the collective flow of various hadrons can be coherently 
explained. 

One can see from the energy dependence of the net- 
baryon rapidity distribution 113911 that the collisions at top 
SPS energy are neither full transparent nor full stopping. 
The suppositions of nuclear stopping power in the two 
models are inappropriate to the nucleus-nucleus collisions 
at top SPS energy. For the description of the rapidity dis- 
tribution for constituent quarks, one can limit the boost in- 
variance into a finite rapidity range in the framework of 
Bjorken model. This modification is often used to analyze 
the longitudinal collectivity in hadronic level ll44ll59ll . The 
rapidity distribution in a isotropic, thermalized fluid ele- 
ment moving with a rapidity 77 is 



X 



1 T I (7) 

cosh(y-i]) cosh-{y-ij) J' ^ 



The rapidity distribution of constituent quarks in the hot 
and dense quark matter is the longitudinal boost-invariant 
superposition of multiple isotropic, thermalized fluid ele- 
ments 



dN,, 
dy 



-(y-T])dT], 



(8) 



rjiiiax is the maximal boot rapidity of fluid elements. The 
average longitudinal collective velocity is taken to be < 
/3l >= tanh(/7,„„.,/2). 

Here, Tj is the temperature of the locally- thermalized 
hot and dense quark matter at hadronization. It is taken 
to he Tf - 170 MeV. m is the constituent mass of quarks 
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when they evolve to the transition point. It is taken to be 
340 MeV for Ught quarks and 500 MeV for strange quarks. 
We have mentioned in above section that the net-quarks, 
still carrying a fraction of initial collision energy, have a 
more complex evolution than hydrodynamic expansion in 
longitudinal axis ll60ll . Therefore, we extract the longitu- 
dinal collective flow from the rapidity distribution of new- 
born quarks. Since most of the data are measured in the 
rapidity range about [-1.5, 1.5], the rapidity spectra of con- 
stituent quarks extracted from experimental data are valid 
only in this region. Using above equations to fit the rapid- 
ity distribution of newborn constituent quarks in Fig. 1, we 
obtain < Pi >- 0.58 for light quarks and < /3l >- 0.65 for 
strange quarks. It is interesting to find that the average lon- 
gitudinal collective velocity of strange quarks is obviously 
greater than that of light quarks. 

For the transverse expansion of the hot and dense quark 
matter, we adopt a blast-wave model proposed by Heniz 
iIsqI] within the boost-invariant scenario. The quarks and 
antiquarks transversely boost with a flow velocity /3r{r) as 
a function of transverse radial position r. /3r{r) is param- 
eterized by the surface velocity /3s: firir) - Ps^", where 
^ = r/R,„ax, and R„ax is the thermal source maximum ra- 
dius (0 < ^ < 1). The transverse momentum distribution 
of constituent quarks in the hot and dense quark matter can 
be equivalently described by a superposition of a set of 
thermalized fluid elements, each boosted with transverse 
rapidity p = tanlr'^Pr 



dN 



2n pjdpT 



r' I pt sinhn\ /rriT coshp\ 

(9) 

Here, Iq and Ki are the modified Bessel functions, mj - 
^Jp^^ + trP- is the transverse mass of the constituent quark. 
The average transverse velocity can be written as 



(fir) 



n + 2' 



Ps. 



(10) 



With fixed parameter n — 0.3, the average transverse ve- 
locity ifS,) is able to characterize the transverse collective 
flow of the hot and dense quark matter Using the above 
equations to fit the transverse momentum distributions of 
the newborn quarks in Fig. 3, we obtain {/?, ) = 0.41 for 
strange quarks and {fir) = 0.36 for light quarks. One can 
see that the (fii) of strange quarks is obviously greater than 
that of light quarks. 

Both longitudinal and transverse results at top SPS en- 
ergy show that the strange constituent quarks get a stronger 
collective flow than the light quarks in the hydrodynamic 
evolution of partonic matter By analyzing the data of 
multi-strange hadrons ll26ll and primary charged hadrons 
the same property is found also at top RHIC energy. 
It suggests that the hot and dense quark matter produced 
at top SPS energy undergoes a similar hydrodynamic evo- 
lution to that at RHIC energies. It is generally believed 
that the decoupled quark a nd g luon plasma (QGP) has been 
created at RHIC energies ||6l|]. This similarity of collective 



flow property in quark level may be regarded as a signal of 
QGP creation at top SPS energy. 



B. The enhanced strangeness 

An interesting phenomenon in high energy heavy ion 
collisions is the enhanced production of strange hadrons, 
which is absent in elementary particle collisions. In rela- 
tivistic heavy ion collisions, enormous amounts of energy 
are deposited in the collision region to create a deconfined 
hot and dense quark matter. The multiple scatterings be- 
tween partons in the hot and dense quark matter will cause 
the large production rate of strangeness by gg — » ss ll27ll . 
The high strangeness of the hot and dense quark matter, af- 
ter hadronization, finally leads to the abundant production 
of the strange hadrons. This phenomenon is regarded as 
a signal of QGP creation. As we know, the enhancement 
of strangeness production at top RHIC energy is quite ob- 
vious [62], and it is generally believed that the QGP has 
been created at RHIC energies. When the collision energy 
drops to the SPS and AGS energies, it is found that the 
strangeness production peaks at about 30A GeV and turns 
into a plateau at higher collision energies ifTill . It is an in- 
dication of the onset of deconfinement. 

In our model, the strangeness of the hot and dense quark 
matter is characterized by the suppression factor As - Ns : 
Nu - Ns : Nj, i.e. the ratio of s quark number to new- 
born u (or d) quark number By fitting the experimental 
data of identified hadrons, we use the model to extract the 
As of hot and dense quark matter at midrapidity in cen- 
tral AA collisions at four energies y/sj^ = 17.3, 62.4, 
130 and 200 GeV, and the results are shown in Table Hill 
The data of midrapidity dN/dy and the calculated results 
with minimum deviations at different collision energies are 
shown also. The statistical uncertainty of As is fixed by 
the twice minimum deviation. The model reproduces the 
hadron yield in reasonably good way, and the chi-square 
fit seems to indicate that with increasing collision energy 
the agreement with data significantly improves. 

One can see that As in such a broad energy range is 
nearly unchanged within statistical uncertainties, exhibit- 
ing an obvious saturation phenomenon. The results of As 
are consistent with the grand canonical limit 0.45) of 
strangeness ifTsll . Using the Bjorken model, one can esti- 
mate that the primordial spatial energy density of the hot 
and dense quark matter produced in collisions at top RHIC 
energy is about 6.0 GeVlfnv' [HI, double of that in Pb-nPb 
collisions at top SPS energy. The difference in primordial 
energy density is large while the final strangeness is nearly 
equal. The hot and dense quark matter created in heavy ion 
collisions is shown to be very close to a perfect fluid iFtiIi . 
It means that the local relaxation time toward to thermal 
equilibrium is much shorter than the macroscopic evolu- 
tion time of the hot and dense quark matter. When the hot 
and dense quark matter evolves to the point of hadroniza- 
tion, the strangeness abundance should be mainly deter- 
mined by the current temperature, irrelevant to the initial 
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TABLE ni: The strange suppression factor /I, and the calculated hadron dN/dy at midr apidity in central AA collisions at different 
energies. The experimental data are taken from Ref. (11 [11 [11 111 111 111 111 Sill 111 S [72]. 



PbPb 17.3 GeV 


Au Au 62.4 GeV 


Au Au 130 GeV 


Au Au 200 GeV 




data 


model 


data 


model 


data 


model 


data 


model 


71+ 


170.1 ±0.7 ±9 


168.3 


212 ± 5.8 ± 14 


211 


276±3±35.9 


268.7 


286.4 ± 24.2 


287.3 


iC 


175.4 ±0.7 ±9 


183.5 


204 ± 7.4 ± 14 


217 


270 ±3.5 ±35.1 


272.4 


281.8 ±22.8 


288.3 


K+ 


29.6 ± 0.3 ± 1.5 


27.3 


33.35 ±2.15 


36.3 


46.7 ± 1.5 ±7 


46.6 


48.9 ±6.3 


48.35 


K- 


16.8 ±0.2 ±0.8 


15.7 


28.16 ± 1.76 


29.9 


40.5 ± 2.3 ± 6 


43.1 


45.7 ± 5.2 


46.48 


P 


29.6 ± 0.9 ± 2.9 


25.9 


27 ± 1.8 ±4.6 


26.17 


19.3 ±0.6 ±3.3 


16.45 


18.4 ±2.6 


17.41 


P 


1.66 ±0.17 ±0.16 


1.53 


11.5 ± 1.5 ±2.9 


11.15 


13.7 ±0.7 ±2.3 


11.63 


13.5 ± 1.8 


13.48 


A 


9.5 ± 0.1 ±1 


13.3 


14.9 ± 0.2 ± 1.49 


13.42 


17.3 ± 1.8 ±2.7 


14.99 


16.7 ± 0.2 ± 1.1 


15.76 


A 


1.24 ±0.03 ±0.13 


1.35 


8.02 ±0.11 ±0.8 


6.77 


12.7 ± 1.8 ±2 


11.4 


12.7 ±0.2 ±0.9 


12.6 




1.44 ±0.1 ±0.15 


1.43 


1.64 ±0.03 ±0.014 


1.63 


2.04 ±0.14 ±0.2 


1.99 


2.17 ±0.06 ±0.19 


2.12 




0.31 ±0.03 ±0.03 


0.26 


0.989 ± 0.057 ± 0.057 


0.96 


1.74 ±0.12 ±0.17 


1.67 


1.83 ±0.05 ±0.2 


1.72 


n + a 




0.17 


0.356 ±0.046 ±0.014 


0.369 


0.56 ±0.11 ±0.05 


0.551 


0.53 ± 0.04 ± 0.04 


0.539 


xVndf 


10.7/7 




6.2/8 




1.6/8 




0.88/8 






0.48 ± 0.09 




0.44 ± 0.02 




0.44 ± 0.04 




0.42 ± 0.025 





energy density and temperature. The same strangeness is 
an indication of the universal hadronization temperature 
for the hot and dense quark matter with low baryon chem- 
ical potential. 

VI. SUMMARY 

In this paper, we have systematically studied the lon- 
gitudinal and transverse production of various hadrons at 
top SPS energy in the scheme of quark combination. Us- 
ing the quark combination model, we firstly calculate the 
yields and rapidity distributions of various hadrons in most 
central Pb+Pb collisions at yfsim = 17.3 GeV. The calcu- 
lated results are in agreement with the experimental data. 
This indicates that the quark combination mechanism is 
applicable in describing the longitudinal hadron produc- 
tion at this collision energy. Secondly the pr distributions 
of various hadrons at top SPS energy are calculated and 
compared with the data. It is found that the light, sin- 
gle and multi-strange hadrons are well reproduced by the 
same quark distributions. It indicates that the hadroniza- 
tion of the hot and dense quark matter is a rapid process. 
The well reproduced baryon/meson ratios in the interme- 



diate pt range at different collision energies are indicative 
of the universality for the quark combination mechanism. 
By fitting the extracted constituent quark distributions at 
hadronization with the hydrodynamic scenario, we further 
obtain the longitudinal and transverse collective flow of the 
hot and dense quark matter produced at top SPS energy. 
It is found that the strange quarks get a stronger collec- 
tive flow than light quarks, which is consistent with that at 
RHIC energies. The strangeness in the hot and dense quark 
matter produced at yfs^ = 17.3, 62.4, 130, 200 GeV are 
extracted. The almost unchanged strangeness may be as- 
sociated with a universal hadronization temperature for the 
hot and dense quark matter with low baryon chemical po- 
tential. 
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